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Abstract. Physical and wind properties of Galactic B supergiants are pre- 
sented based upon non-LTE line blanketed model atmospheres, including Sher 25 
toward the NGC 3603 cluster. We compare Ha derived wind densities with re- 
cent results for SMC B supergiants and generally confirm theoretical expecta- 
tions for stronger winds amongst Galactic supergiants. Mid B supergiant winds 
are substantially weaker than predictions from current radiatively driven wind 
theory, a problem which is exacerbated if winds are already clumped in the Ha 
line forming region. We find that the so-called 'Instability jump' at Bl (T g 
~ 21kK) from Lamers et al. is rather a more gradual downward trend. CNO 
elemental abundances, including Sher 25, reveal partially processed material at 
their surfaces. In general, these arc in good agreement with evolutionary pre- 
dictions for blue supergiants evolving redward accounting for rotational mixing. 
A few cases, including HD 152236 (C 1 Sco), exhibit strongly processed material 
which is more typical of Luminous Blue Variables. Our derived photospheric 
[N/O] ratio for Sher 25 agrees with that for its ring nebula, although a higher 
degree of CNO processing would be expected if the nebula originated during a 
red supergiant phase, as is suspected for the ring nebula ejected by the B su- 
pergiant progenitor of SN 1987A, Sk-69° 202. Sher 25 has an inferred age of 
~5Myr in contrast with ~2Myr for HD 97950, the ionizing cluster of NGC 3603, 
so it may be a foreground object or close binary evolution may be responsible 
for its unusual location in the H-R diagram. 



1. Introduction 

O stars dominate the energetics of young starbursts since their ionizing output 
is a very steep function of effective temperature (Leitherer et al. 1992). As 
such, numerous intensive studies of O stars in the Milky Way and Magellanic 
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Clouds have been undertaken, with regard to determining their fundamental 
parameters and the empirical dependence of their wind strengths on metallicity 
(e.g. Mokiem et al. 2006ab). Unfortunately, O dwarfs are visually faint in 
external galaxies due to their large bolometric corrections and small radii, with 
typically My = -5.0 mag. O star abundances of CNO elements are difficult to 
derive (Crowther et al. 2002), yet these are sensitive to early rotational mixing 
(e.g. Meynet & Maeder 2000). 

In contrast, B supergiants have received rather less attention since they are 
less relevant for energetics of young starbursts. Nevertheless, they are visually 
bright, typically My = —7.0 mag, possess larger radii and so are easy to observe 
individually in nearby galaxies, permitting robust tests of the metallicity de- 
pendence of radiatively driven wind theory, which is predicted to increase from 
early to mid B supergiants (Vink et al. 2000). In addition, CNO abundances 
are readily obtained from optical spectroscopy, for comparison to evolutionary 
model predictions. Notably, Trundle et al. (2004) and Trundle & Lennon (2005) 
have studied a large sample of SMC B supergiants using the line blanketed ver- 
sion of FASTWIND (Puis et al. 2005). To date, the most extensive study 
of Galactic B supergiants is by Kudritzki et al. (1999), who adopted the T e g 
scale from the plane-parallel unblanketed study of McErlean et al. (1999), plus 
mass- loss rates from the unblanketed version of FASTWIND. 

Our present study of Galactic B supergiants employs CMFGEN (e.g. Hillier 
et al. 2003) which also incorporates line blanketing and spherical geometry, and 
so permits a direct comparison with recent SMC results, plus CNO abundance 
determinations. Do the abundances of normal B supergiants match those of 
stars evolving towards red supergiants (RSG) or returning from RSG following 
convective dredge-up? Comparisons with Luminous Blue Variables (LBVs) are 
also made, since they exhibit spectral characteristics of early B supergiants at 
visual minimum and A supergiants at visual maximum (Crowther 1997). No- 
tably, C 1 Sco (HD 152236) is an early B hypergiant and lies above the empirical 
Humphreys-Davidson limit (Humphreys h Davidson 1979). We also include 
new stellar and nebular studies of Sher 25 (Hendry et al. 2006), an apparently 
normal B supergiant in the giant HII region NGC 3603 yet possessing an ejecta 
nebula reminiscent of that associated with SN 1987A (Brandner et al. 1997b). 

Finally, we consider the empirical bistability jump amongst early B super- 
giants, which was originally developed for the LBV P Cygni by Pauldrach & Puis 
(1990), in the sense that slight variations in its effective temperature resulted 
in hydrogen being ionized or neutral in its outer stellar wind (Najarro et al. 
1997). Lamers et al. (1995) studied wind velocities of OB stars, and established 
a discontinuity in the ratio of to v csc for B0. 5 and Bl supergiants, which is 
assumed in current stellar wind models of Vink et al. (2000). 

2. Observations and analysis 

Full details of observations and method of analysis are presented by Crowther et 
al. (2006). In summary we have observed 26 Galactic early- type la supergiants 
- spanning the spectral range 09.5 to B3 - at intermediate dispersion in the blue 
and red using the 1.0m JHK, 2.5m INT or 4.2m WHT for the northern sample 
or at the 1.5m CTIO for the southern sample, with the exception of Sher 25 
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Figure 1. Temperature scale of B supergiants from this Milky Way study 
(CMFGEN, circles) plus Trundle et al. (2004) and Trundle & Lennon (2005) 
for the SMC (FASTWIND, triangles) versus the McErlean et al. (1998, solid 
line) calibration. 



for which 3.9m AAT echelle spectroscopy was obtained. Cluster or association 
distance estimates are available for 75% of the sample, whilst average subtype 
absolute magnitudes are taken from Magellanic Cloud stars for the remainder. 
For the case of Sher 25, we adopt a distance of 7.9 kpc (Russeil 2003) towards 
NGC 3603 from which we obtain My = —7.4 mag. 

For the present study we have used CMFGEN (see Hillier et al. 2003) 
which solves the radiative transfer equation in the co-moving frame, under the 
constraints of radiative and statistical equilibrium. For the supersonic part, the 
velocity is parameterized with a classical (3 law, which is connected to a hy- 
drostatic density structure at depth using a H-He TLUSTY model. A depth 
independent Doppler profile was assumed for the atmospheric structure calcu- 
lation in the co-moving frame, whilst a radially dependent turbulence was used 
to calculate the emergent spectrum in the observer's frame, in which incoherent 
electron scattering and Stark broadening for H, He were included. 

Stellar temperatures were derived from suitable lines of Hel-II (O9.5-B0), 
Silll-IV (B0.5-B2), Sill-III (B2.5-B3) such that radii and luminosities followed 
from absolute magnitudes, with Ha observations providing M and the veloc- 
ity exponent (3. Terminal velocities and rotational velocities were taken from 
Howarth et al. (1997). H/He abundances are difficult to accurately derive in 
B supergiants, such that He/H=0.2 by number is adopted throughout. CNO 
abundances were varied to reproduce the relevant optical absorption lines, and 
compared with Solar abundances (Asplund et al. 2004). 
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Figure 2. Observed H-R diagram for our sample of B supergiants, coded 
by their degree of N/O enrichment, together with tracks from Solar metal- 
licity, rotating (t>i n i t =300 km sec -1 ) evolutionary models (dotted lines) from 
Meynet & Maeder (2000) for which predicted N/O enrichments are indicated 
for redward (blueward) evolution at T e ff = 22kK. The location of Sher 25 and 
HD 152236 (C 1 Sco) are indicated. 



3. Physical and Wind Properties 
3.1. Temperature scale 

The inclusion of line blanketing and spherical geometry into O star models has 
resulted in a major 15±5% downward revision in their T e g- spectral type calibra- 
tion (Crowther et al. 2002; Martins et al. 2005). For B supergiants, we present 
our results together with those of Trundle et al. (2004) and Trundle & Lennon 
(2005) for SMC stars studied with FASTWIND (Puis et al. 2005) in Fig. [TJ 
together with the results from unblanketed plane parallel models of McErlean 
et al. (1999). As is apparent, the revision in B supergiant temperature scale is 
modest, with typical revisions of —1 to — 2kK. Note that results for line blan- 
keted versions of CMFGEN and TLUSTY are consistent to within ±0.5kK for 
several stars in common between the present study and Urbaneja (2004). 

The location of our stars in the H-R diagram is presented in Fig. |5J such 
that the majority are consistent with initial masses in the range 20-40 Mq, 
whilst C 1 Sco is rather more massive. Spectroscopically derived masses span a 
wide range, from 8 up to 76 Mq, based upon surface gravities reliable to logg 
± 0.15 dex. In most cases spectroscopic masses are lower than those estimated 
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from evolutionary models, indicating another example of a 'mass discrepancy' 
in common with Trundle & Lennon (2005). 
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Figure 3. Comparison between the wind momenta of Galactic BO. 7-3 su- 
pergiants ( T c ff <23kK) and studies of SMC counterparts Trundle et al. 2004; 
Trundle & Lennon 2005), assuming homogeneous Ha mass- loss rates, plus 
radiatively driven wind theory predictions by Vink et al. (2000, 2001) for the 
Milky Way (solid) and SMC (dotted). Sher 25 is indicated, revealing wind 
properties typical of Galactic B supergiants. 



3.2. Mass-loss rates 

Ha mass-loss rates compare well with mid-IR excess methods (e.g. Barlow 
& Cohen 1977) and line blanketed FASTWIND models (Urbaneja 2004), but 
exceed those derived by Kudritzki et al. (1999) by a factor of three, who used an 
unblanketed version of FASTWIND for B2-3 supergiants. Spectroscopic studies 
in which clumped winds are included tend to produce superior spectral fits to 
far-UV lines (Evans et al. 2004b), suggesting comparable clumping factors for 
the inner (Ha line) and outer (mid-IR continuum) forming regions. 

In general, wind momenta of Galactic BO-0.5 subtypes lie close to the Vink 
et al. (2000) prediction. If winds of early B supergiants are clumped in the Ha 
line forming region, we would need to shift the observed wind momenta to lower 
values. Fig. El presents wind momenta for Galactic B0. 7-3 supergiants (T e ff (< 
23kK) together with corresponding SMC results from Trundle et al. (2004) and 
Trundle & Lennon (2005), plus theoretical predictions from Vink et al. (2000, 
2001). A clear separation between Galactic and SMC supergiants is appar- 
ent, with approximately the offset predicted by Vink et al. (2001). However, 
measured values lie typically ~0.5 dex below the Vink et al. calibration. As dis- 
cussed above, if Ha derived mass-loss rates of OB stars need to be corrected for 
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clumped winds, as suggested by recent observational evidence, including Evans 
et al. (2004b) for B supergiants, the difference between theory and observation 
would be further exacerbated. 
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Figure 4. The ratio Woo/w eS c as a function of effective temperature for Galac- 
tic B supergiants from the present study (filled circles), together with results 
from Evans et al. (2004a, open triangles) for SMC B supergiants. Error bars 
relate to uncertainties in temperature (±1 kK), gravity (±0.15 to 0.2 dex) 
and absolute magnitude (±0.3 to 1 mag). 



3.3. Bistability jump 

In Fig. 0] we present our current results, together with values from Evans et al. 
(2004a), based on measured HST/STIS wind velocities, plus model atmosphere 
results from Evans et al. (2004b), Trundle et al. (2004), Trundle & Lennon 
(2005). In contrast with the discontinuity identified by Lamers et al. (1995), we 
identify a gradual downward trend of Voq / v csc with temperature, albeit with 
a large scatter. For T c g> 24kK (approximately B0. 5 la and earlier), / v csc 
~ 3.4, for 20kK< T e g <24kK (approximately B0. 7-1 la) / v csc ~ 2.5, and for 
T c ff <20kK (approximately B1.5Ia and later) Voo / v esc ~ 1.9. This reveals that 
the Bl 'jump' is misleading in the context of normal B supergiants. Prinja &: 
Massa (1998) came to similar conclusions based on a larger B supergiant sample, 
albeit with an adopted subtype-temperature calibration. 

Empirical values of the Lamers et al. (1995) 'bistability jump' were adopted 
by Vink et al. (2000) in their radiatively driven wind calculations. These pre- 
dictions have subsequently been used in evolutionary models (e.g. Meynet & 
Maeder 2000). Consequently, our re-determination of physical parameters and 
wind properties of early B supergiants has potential consequences for evolution- 
ary and spectral synthesis calculations. 
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3.4. Photospheric CNO abundances 

Fig. |S] provides a summary of our derived CNO abundances. With respect to 
the current Solar values (Asplund et al. 2004), carbon is depleted by 0.4+0.4 
dex, oxygen is mildly depleted by up to 0.5 dex, whilst nitrogen is enriched 
by 0.6+0.4 dex. Mean [N/C] and [N/O] ratios indicate evidence for partial 
CNO processing within the photospheres of B supergiants. Our results support 
previous abundance estimates from the literature, which were typically based 
upon differential abundance analyses from plane-parallel LTE or non-LTE model 
atmospheres. 

On average, the observed degree of CNO processing agrees well with recent 
evolutionary models for 2O-4OM , for which [N/C] = +0.8 to +1.3 and [N/O] 
= +0.7 to +1.2 were predicted via rotational mixing (wmit = 300 km sec -1 ) at 
~22kK on the redward evolutionary track (Meynet & Maeder 2001). Spectro- 
scopically, all our sample of supergiants are morphologically normal, except for 
k Cas (HD 2905, BC0.7Ia) which indeed possesses the lowest N enrichment in 
our sample with [N/O] = +0.1 dex. Conversely, C 1 Sco (HD 152236, B1.5 Ia+), 
one of two hypergiants in our sample, possesses the highest N enrichment of 
[N/O] = +1.3, that is more typical of nebulae associated with Luminous Blue 
Variables, e.g. [N/O] = +1.3 to +1.9 dex for AG Car and R127 (Lamers et al. 
2001). 



4. Sher 25 in NGC 3603 

Sher 25 (B1.5 lab, Moffat 1983) is associated with a ring nebula and apparent 
bipolar outflows (Brandner et al. 1987a). The ring nebula is reminiscent of the 
inner nebula associated with SN 1987A, which is presumed to have been ejected 
from the Sk-69° 202 (B3 I, Walborn et al. 1989) progenitor. Both nebulae 
share physical sizes of ~0.5 pc, with dynamical ages of several thousand years 
and apparent nitrogen enrichment (Panagia et al. 1996; Brandner et al. 1997b). 
Hendry et al. (2006) present studies of the ring nebula and bipolar lobes from 
Sher 25, plus that of the giant HII region NGC 3603 toward which Sher 25 is 
located. 

The stellar and wind properties of Sher 25 are typical of other Galactic 
early B supergiants in our sample (recall Figs. [2H21)> as derived from opti- 
cal AAT/UCLES spectroscopy, i.e. T cS = 21kK, log L* / L e = 5.68, M = 
1.2xl0" 6 M Q yx- 1 Voo ^ 750 km sec 1 using a distance of 7.9 kpc (Russeil 2003). 
The surface elemental abundances of Sher 25 derived here, [N/C]=+1.3 and 
[N/0]=+1.2, are compared with other B supergiants in Fig. [HJ Table ^provides 
a comparison with previous photospheric derivations (Smartt et al. 2002) plus 
ring nebula and bipolar lobe abundances in Tabled Note that our derived [N/O] 
ratio agrees perfectly with the mean nebular ratio from Hendry et al. (2006), 
whilst the giant HII region is close to Solar composition (see also Esteban et al. 
2005). 

If we compare the location of Sher 25 in the H-R diagram with evolutionary 
models of Meynet & Maeder (2000) it lies close to an initial 4OM rotating with 
v mit = 300 km sec -1 (recall Fig. [2J). If we adjust the predictions from Meynet 
& Maeder to the Solar abundances from Asplund et al. (2004) one expects 
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Figure 5. [N/C] (left) and [N/O] (right) abundances for Milky Way B su- 
pcrgiants from the present CMFGEN study (relative to Asplund ct al. 2004), 
typically revealing partially CNO processed material, in which Sher 25 is 
indicated in grey. Although the relative abundances presented here are an- 
ticipated to be robust, there remains considerable absolute uncertainty, since 
results for a subset of the present sample using FASTWIND differ by up to 
0.5 dex (Urbaneja 2004). 



a nitrogen enhancement of +1.0 dex and oxygen depletion of -0.2 dex on the 
redward track after 5.1 Myr, i.e. [N/O] = +1.2 dex, in perfect agreement with 
the mean nebular and stellar CMFGEN results. For the post-RSG stage after 
5.2 Myr, [N/O] = +1.8 dex is predicted. Consequently, Sher 25 is consistent 
with a single star evolving towards the RSG phase. 

However, amongst the most successful models to explain SN 1987A during 
its B supergiant phase are those of Podsiadlowski (1992), involving either ac- 
cretion from, or merger with, a binary companion during the common envelope 
evolution in the RSG stage of Sk-69° 202. Assuming membership of NGC 3603, 
Sher 25 is more luminous than Sk-69° 202 - for which log L*j £©~ 5.1. If 
Sher 25 has undergone a similar evolution, which may be suggested from its 
recent nebular ejection, this star would not be expected to possess CNO abun- 
dances typical of normal B supergiants evolving towards the RSG phase. The 
dynamical age of the ring nebula and bipolar lobes is ~8700 yr (Brandner et al. 
1997b, adjusted for a distance of 7.9 kpc). 

One final puzzle regarding Sher 25 is that it lies only 20 arcsec (1 pc) 
in projection from the central ionizing cluster, HD 97950, in NGC 3603, yet 
its inferred age of ~5Myr greatly exceeds that of HD 97950 which is ~2Myr 
(Crowther & Dessart 1998). HD 97950 closely resembles R136 in the LMC 
(Walborn 1973; Moffat 1983), yet unlike R136, it does not possess an extensive 
halo of earlier phases of massive star formation. Perhaps Sher 25 is a foreground 
object, which may plausibly be inferred from the presence of other stars which 
are older than the HD 97950 cluster to one side (Brandner et al. 1997a). In 
this case the stellar luminosity and the size of the ring nebula are smaller than 
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adopted here, although we note its physical and wind properties are typical of 
other Galactic early B supergiants. Alternatively, close binary evolution may be 
responsible for its unusual position in the H-R diagram. 

Table 1. Photosphere abundances of Sher 25, together with nebular abun- 
dances for the ejecta nebula and NGC 3603 HII region, with reference to Solar 
abundances of log(0/H) + 12=8.66 and log(N/H)+12=7.78 from Asplund et 
al. (2004). 
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log N/H 


[N/O] 
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Star (TLUSTY) 
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Ring nebula 
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Hendry et al. 2006 
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